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Abstract. We develop a statistical theory that describes quantum-mechanical 
scattering of a particle by a cavity when the geometry is such that the classical 
dynamics is chaotic. This picture is relevant to a variety of physical systems, rang- 
ing from atomic nuclei to mesoscopic systems and microwave cavities; the main 
application to be discussed in this contribution is to electronic transport through 
mesoscopic ballistic structures or quantum dots. The theory describes the regime 
in which there are two distinct time scales, associated with a prompt and an equili- 
brated response, and is cast in terms of the matrix of scattering amplitudes S. We 
construct the ensemble of S matrices using a maximum-entropy approach which 
incorporates the requirements of flux conservation, causality and ergodicity, and 
the system-specific average of S which quantifies the effect of prompt processes. 
The resulting ensemble, known as Poisson's kernel, is meant to describe those situ- 
ations in which any other information is irrelevant. The results of this formulation 
have been compared with the numerical solution of the Schrodinger equation for 
cavities in which the assumptions of the theory hold. The model has a remarkable 
predictive power: it describes statistical properties of the quantum conductance of 
quantum dots, like its average, its fiuctuations, and its full distribution in several 
cases. We also discuss situations that have been found recently, in which the notion 
of stationarity and ergodicity is not fulfilled, and yet Poisson's kernel gives a good 
description of the data. At the present moment we are unable to give an explana- 
tion of this fact. 

Keywords: Quantum Chaotic Scattering, Statistical S matrix, Mesoscopic Physics, 
Information Theory, Maximum Entropy. 
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1 Introduction 

The problem of eoherent multiple scattering of waves has long been of great 
interest in physics and, in particular, in optics. There are a great many 
wave-scattering problems, appearing in various fields of physics, where the 
interference pattern due to multiple scattering is so complex, that a change in 
some external parameter changes it completely and, as a consequence, only 
a statistical treatment is feasible and, perhaps, meaningful [l|2j . 

In the problems to be discussed here, complexity in wave scattering de- 
rives from the chaotic nature of the underlying classical dynamics. Our 
discussion will find applications to problems like: i) electronic transport in 
microstructures called ballistic quantum dots, and ii) transport of electro- 
magnetic waves, or other classical waves (like elastic waves), through cavities 
with a chaotic classical dynamics. In particular, we shall study the statistical 
fluctuations of transmission and reflection of waves, which are of considerable 
interest in mesoscopic physics. 

The ideas involved in our discussion have a great generality. Let us re- 
call that, historically, nuclear physics -a complicated many-body problem- 
has offered very good examples of complex quantum-mechanical scattering. 
Here, the typical dimensions are on the scale of the fm (1 fm = 10^^'^cm). 
The statistical theory of nuclear reactions has been of great interest for many 
years, in those cases where, due to the presence of many resonances, the cross 
section is so complicated as a function of energy that its detailed structure 
is of little interest and a statistical treatments is then called for. Most re- 
markably, one finds similar statistical properties in the quantum-mechanical 
scattering of "simple" one-particle systems -like a particle scattering from a 
cavity- whose classical dynamics is chaotic. The typical dimensions of the 
ballistic quantum dots that were mentioned above is 1 fim, while those of mi- 
crowave cavities is of the order of 0.5 meters. Thus the size of these systems 
spans ~ 14 orders of magnitude! 

The purpose of this contribution is to review past and recent work in 
which various ideas that were originally developed in the framework of nuclear 
physics, like the nuclear optical model and the statistical theory of nuclear re- 
actions [3] , have been used to give a unified treatment of quantum- mechanical 
scattering in simple one-particle systems in which the corresponding classical 
dynamics is chaotic, and of microwave propagation through similar cavities 
(see Ref. [1 and references contained therein). 

The most remarkable feature that we shall encounter is the statistical 
regularity of the results, in the sense that they will be expressible in terms 
of a few relevant physical parameters, the remaining details being just "scaf- 
foldings". This feature will be captured within a framework that we shall 
call the "maximum-entropy approach" . 

In the next section we present some general physical ideas related to 
scattering of microwaves by cavities and of electrons by ballistic quantum 
dots. The formulation of the scattering problem is given in Sec. [3l where we 
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employ, to be specific, the language of electron scattering. The notion of do- 
ing statistics with the scattering matrix and the resulting maximum-entropy 
model are reviewed in Sees. |4]and[5l respectively. In Sec. [6] we compare the 
predictions of our model with a number of computer simulations. We also 
present a critical discussion of the results, based on some recent findings on 
the the validity conditions of our statistical model. The conclusions of this 
presentation are the subject of Sec. [71 

2 Microwave cavities and ballistic quantum dots 

Consider the system shown schematically in Fig. [1] It consists of a cavity 
connected to the external world by two waveguides, ideally of infinite length. 




Fig. 1. The 2D cavity referred to in the text. The cavity is connected to the 
outside via two waveguides, each supporting A'^ running modes. The arrows inside 
the waveguides indicate incoming or outgoing waves. In waveguide 1 — 1,2 the 
wave amplitudes are indicated as an\ fen\ respectively, where n = 1, • • ■, A''. 

We may think of performing a microwave experiment with such a device 
[3]. A wave is sent in from one of the waveguides, the result being a reflected 
wave in the same waveguide and a transmitted wave in the other. One 
given point inside the system receives the contributions from many multiply 
scattered waves that have bounced from the inner surface of the cavity, giving 
rise to the complex interference pattern mentioned in the Introduction. An 
alternative interpretation of the same situation is that the outgoing wave has 
suffered the effect of many resonances that are present inside the cavity. The 
net result is an interference pattern which shows an appreciable sensitivity to 
changes in external parameters: for example, a small variation in the shape 
of the cavity may change the pattern drastically. 

In experiments on electronic transport [5^ performed with ballistic mi- 
crostructures, or quantum dots, one finds, for sufficiently low temperatures 
and for spatial dimensions on the order of l/xm or less, that the phase coher- 
ence length l^ exceeds the system dimensions; thus the phase of the single- 
electron wave function -in an independent-electron approximation- remains 
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coherent across the system of interest. Under these conditions, these systems 
are cahed mesoscopic |6|7|8j . The elastic mean free path l^i also exceeds the 
system dimensions: impurity scattering can thus be neglected, so that only 
scattering from the boundaries of the system is important. In these sys- 
tems the dot acts as a resonant cavity and the leads as electron waveguides. 
Experimentally, an electric current is established through the leads that con- 
nect the cavity to the outside, the potential difference is measured and the 
conductance G is then extracted. Assume that the microstructure is placed 
between two reservoirs (at different chemical potentials) shaped as expand- 
ing horns with negligible reflection back to the microstructure. In a picture 
that takes into account the electron-electron interaction in the random-phase 
approximation (RPA), the two-terminal conductance at zero temperature - 
understood as the ratio of the current to the potential difference between the 
two reservoirs- is then given by Landauer's formula (101911] 

G = 2^T, (1) 

where the total transmission T is given by 

T = tr (ttt) , (2) 

t being the matrix of transmission amplitudes tab associated with the resulting 
self-consistent single-electron problem at the Fermi energy cp, with a and b 
denoting final and initial channels. The factor of 2 in Eq. ^ is due to 
the two-fold spin degeneracy in the absence of spin-orbit scattering that we 
shall assume here. In this "scattering approach" to electronic transport, 
pioneered by Landauer [10], the problem is thus reduced to understanding 
the quantum-mechanical single-electron scattering by the cavity under study. 

When the system is immersed in an external magnetic field B and the 
latter is varied, or the Fermi energy e^? or the shape of the cavity are varied, 
the relative phase of the various partial waves changes and so do the interfer- 
ence pattern and the conductance. This sensitivity of G to small changes in 
parameters through quantum interference is called conductance fluctuations. 

The experiments on quantum dots given in Ref. [5] report cavities in the 
shape of a stadium, for which the single-electron classical dynamics would 
be chaotic, as well as experimental ensembles of shapes. The average of the 
conductance, its fluctuations and its full distribution were obtained over such 
ensembles. 

In what follows we shall phrase the problem in the language of electronic 
scattering; spin will be ignored, so that we shall deal with a scalar wave 
equation. However, our formalism is applicable to scattering of microwaves, 
or other classical waves, in situations where the scalar wave equation is a 
reasonable approximation. 
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3 The scattering problem 



In our application to mesoscopic physics we consider a system of noninteract- 
ing "spinless" electrons and study the scattering of an electron at the Fermi 
energy ei;- by a 2D microstructure, connected to the outside by two leads of 
width W (see Fig. [I]). 

Very generally, a quantum-mechanical scattering problem is described by 
its scattering matrix, or S matrix, that we now define. The confinement of 
the electron in the transverse direction in the leads gives rise to the so called 
transverse modes, or channels. If the incident Fermi momentum kp satisfies 



N < kpW/T: < 1, 



(3) 



there are N transmitting or running modes (open channels) in each of the 
two leads. The wavefunction in lead / {I — 1,2 for the left and right leads, 
respectively) is written as the A^-dimensional vector 



(4) 



the n-th component {n — 1, . . . , N) being a linear combination of plane waves 
with constant flux (given by h), i.e.. 



-iknXi 



+ h 



(27rrfc„/m)i/2 (27rr/fc„/m)i/2 



1, 



(5) 



m being the electron mass. In ([5]), kn is the "longitudinal" momentum in 
channel n, given by 

(6) 



W 



— kp. 



The 2iV-dimensional S matrix relates the incoming to the outgoing ampli- 
tudes as 



6(1) 

6(2) 



= s 



ad) 
a(2) 



(7) 



where ad), a(2), 6(1), 6(2) are iV-dimensional vectors. The matrix S has the 
structure 

'r 
t r' 



(8) 



where r, t are the NxN reflection and transmission matrices for particles from 
the left and r', t' for those from the right: this thus defines the transmission 
matrix t used earlier in Eq. . 

The requirement of flux conservation (FC) implies unitarity of the S ma- 
trix [T], i.e., 

SS^ = /. (9) 

This is the only requirement in the absence of other symmetries. This is the 
so-called unitary case, also denoted by /3 = 2 in the literature on random- 
matrix theory. For a time-reversal-invariant (TRI) problem (as is the case 
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in the absence of a magnetic field) and no spin, the S matrix, besides being 
unitary, is symmetric [1]: 

S = S^. (10) 

This is the orthogonal case, also denoted by /3 = 1. The symplectic case {[3 — 
4) arises in the presence of half-integral spin and time-reversal invariance, 
and in the absence of other symmetries, but will not be considered here. 

4 Doing statistics with the scattering matrix 

In numerical simulations of quantum scattering by 2D cavities with a chaotic 
classical dynamics one finds that the S matrix and hence the various trans- 
mission and reflection coefficients fluctuate considerably as a function of the 
incident energy E (or incident momentum fc), because of the resonances oc- 
curring inside the cavity [1] . These resonances are moderately overlapping for 
just one open channel and become more overlapping as more channels open 
up. An example is illustrated in Fig. [51 which shows the total transmission 




kW/TT 

Fig. 2. Total transmission coefficient of Eq. ((2]| as a function of the incident mo- 
mentum, obtained by solving numerically the Schrodinger equation for the system 
indicated in the upper left corner (from Ref. [H]), whose classical dynamics has a 
chaotic nature. 

coefficient T{E) of Eq. ([2]) [proportional to the conductance, by Landauer's 
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formula of Eq. ([ij], obtained by solving numerically the Schrodinger equa- 
tion for a cavity whose underlying classical dynamics is chaotic. Since we 
are not interested in the detailed structure of T{E), it is clear that we are 
led to a statistical description of the problem, just as was explained in the 
Introduction. 

As E changes, S{E) wanders on the unitarity sphere, restricted by the 
condition of symmetry if TRI applies. The question we shall address is: what 
fraction of the "time" (energy) do we find S inside the "volume element" 
dii{S)7 [dfj,{S) is the invariant measure for S matrices, to be defined below.] 
We call dP{S) this fraction and write it as 

dP{s)^p{S)d^l{s), (11) 

where p{S) will be called the "probability density" . Once we phrase the prob- 
lem in this language, we realize that we are dealing with a random S-matrix 
theory. In the next section we shall introduce a model for the probability 
density p{S). 

The concept of invariant measure dfi{S) for S matrices mentioned above is 
an important one and we now briefly describe it. By definition, the invariant 
measure for a given symmetry class does not change under an automorphism 
of that class of matrices onto itself |12|13j . i.e., 

d^l^>^\s) = d^l^>^\s'). (12) 

For unitary and symmetric 5* matrices we have 

S' = UoSU^, (13) 

Uq being an arbitrary, but fixed, unitary matrix. Eq. represents an 

automorphism of the set of unitary symmetric matrices onto itself and de- 
fines the invariant measure for the orthogonal (/? = 1) symmetry class of S 
matrices. We recall that this class is relevant for systems when we ignore 
spin and in the presence of TRI. When TRI is broken, we deal with unitary, 
not necessarily symmetric, S matrices. The automorphism is then 

S' - Uf^SVo, (14) 

Uq and Vq being now arbitrary fixed unitary matrices. For this unitary, 
or /3 = 2, symmetry class of S matrices the resulting measure is the well 
known Haar measure of the unitary group and its uniqueness is well known. 
Uniqueness for the /3 = 1 class (and /3 = 4, not discussed here) was shown 
in Ref. [12]. The invariant measure, Eq. (fT2|) . defines the so-called Circular 
(Orthogonal, Unitary) Ensemble (COE, CUE) of S matrices, for /3 = 1,2, 
respectively. 

5 The maximum-entropy model 

We shall assume that in the scattering process by our cavities two distinct 
time scales occur |l|3|14ll5|16|17j : i) a prompt response arising from short 
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paths; ii) a delayed response arising from very long paths. HistoricaUy, this 
assumption was first made in the nuclear theory of the compound nucleus, 
where the prompt response is associated with the so called direct processes, 
and the delayed one with an equilibrated response arising from the compound- 
nucleus formation. 

The prompt response is described mathematically in terms of S, the av- 
erage of the actual scattering matrix S{E) over an energy interval around a 
given energy E. Following the jargon of nuclear physics, we shall refer to S 
as the optical S matrix. The resulting averaged, or optical, amplitudes vary 
much more slowly with energy than the original ones. 

As explained below, the statistical distribution of the scattering matrix 
iS'(ii^) in our chaotic-scattering problem is constructed through a maximum- 
entropy "ansatz" , assuming that it depends parametrically solely on the op- 
tical matrix S, the rest of details being a mere "scaffolding" . 

As in the field of statistical mechanics, it is convenient to think of an 
ensemble of macroscopically identical cavities, described by an ensemble of 
S matrices (see also Refs. |18I19| ). In statistical mechanics, time averages 
are very difficult to construct and hence are replaced by ensemble averages 
using the notion of ergodicity. In a similar vein, in the present context we 
idealize S{E), for all real E, as a stationary random-matrix function of E 
satisfying the condition of ergodicity |1I17| . so that we may study energy 
averages in terms of ensemble averages. For instance, the optical matrix S 
will be calculated as an ensemble average, i.e., S — (S), which will also be 
referred to as the optical S matrix. 

We shall assume E to be far from thresholds: locally, S{E) is then a 
meromorphic function which is analytic in the upper half of the complex- 
energy plane and has resonance poles in the lower half plane. 

From "analyticity- ergodicity" one can show the "reproducing" property: 
given an "analytic function" of 5, i.e., 

= Caibf-apbp Sa^bi ■ ■ ■ Sapbp, (15) 

a\b\---apbp 

we must have: 

{f{S)) EE J fiS) p^s) (S) df,{S) = fiiS)). (16) 

This is the mathematical expression of the physical notion of analyticity- 
ergodicity. Thus p(^s) (5*) must be a "reproducing kernel". Notice that only 
the optical matrix (S) enters the definition. The reproducing property (|16p 
and reality of the distribution do not fix the ensemble uniquely. However, 
among the real reproducing ensembles, Poisson's kernel jl3| , i.e.. 



p\s]iS) I det(/-5(5t)) |2W/3+2-/3 ' ^^^^ 



,i[det(/-(5)(5t))](2^/'+2-'3)/2 
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(recall that 2N is the dimensionality of the S matrices, N being the number 
of open channels in each lead) is special, because its information entropy 



is greater than or equal to that of any other probability density satisfying 
the reproducibility requirement for the same optical (5). We could describe 
this result qualitatively saying that, for Poisson's kernel, S is as random as 
possible, consistent with (S) and the reproducing property. 

As for its information-theoretic content, Poisson's kernel of Eq. (|17p de- 
scribes a system with: I.- the general properties associated with i) unitarity of 
the S matrix (flux conservation), ii) analyticity of S{E) implied by causality, 
iii) presence or absence of time-reversal invariance (and spin-rotation sym- 
metry when spin is taken into account) which determines the universality 
class [orthogonal {(3 = 1), unitary (/3 = 2) or symplectic (/? = 4)], and II.- 
the system-specific properties -parametrized by the ensemble average (S)-, 
which describe the presence of short-time processes. System-specific details 
other than the optical (S) are assumed to be irrelevant. 

6 Comparison of the information-theoretic model with 
numerical simulations 

A number of computer simulations have been performed, in which the Schrodinger 
equation for 2D structures was solved numerically and the results were com- 
pared with our theoretical predictions. In Refs. |1I15I16I18] ). statistical 
properties of the quantum conductance of cavities were studied, like its av- 
erage, its fluctuations and its full distribution, both in the absence and in 
the presence of a prompt response. Here we concentrate on the conductance 
distribution, mainly in the presence of prompt processes. 

Fig. [3] shows the probability distribution of the conductance [which is 
proportional to the total transmission coefficient, according to Landauer's 
formula, Eq. ([1])] for a quantum chaotic billiard with one open channel in 
each lead {N — 1), embedded in a magnetic field: the relevant universality 
class is thus the unitary one (/3 — 2) (reported from Refs. |l|15|16j ). The 
numerical conductance distribution was obtained collecting statistics by: 

i) sampling in an energy window AE larger than the energy correlation 
length, but smaller than the interval over which the prompt response changes, 
so that (S) is approximately constant across the window AE. Typically, 200 
energies were used inside a window for which kW/n G [1.6, 1.8] (where W is 
the width of the lead) . That such energy scales can actually be found in the 
physical system under study is interpreted as giving evidence of two rather 
widely separated time scales in the problem. 

ii) sampling over 10 slightly different structures, obtained by changing the 
height or angle of the convex "stopper" , used mainly to increase the statistics. 




(18) 
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Fig. 3. The probability distribution of the total transmission coefficient T, or con- 
ductance, obtained by numerically integrating the Schrodinger equation for the 
billiards shown on the left side of the figure for A'^ = 1. The results are indi- 
cated with squares, that include the statistical error bars. The curves are obtained 
from Poisson's kernel, Eq. (|17p . with (S) extracted from the numerical data. The 
agreement is good. (From Ref. [15].) 



Since we are for the most part averaging over energy, we rely on er- 
godicity to compare the numerical distributions with the ensemble averages 
of the theoretical maximum-entropy model. The optical 5* matrix was ex- 
tracted directly from the numerical data and used as (5) in Eq. pT)l for 
Poisson's kernel; in this sense the theoretical curves shown in the figure rep- 
resent parameter-free predictions. 

Several physical situations were considered in order to vary the amount 
of direct processes, quantified through (S). The upper panels in Fig. [3] 
correspond to structures with leads attached outside the cavity, whereas the 
lower panels show the conductance distribution when the attached leads are 
extended into the cavity; in the latter case, the presence of direct transmission 
is promoted. In some instances, a tunnel barrier with transmission coefficient 
1/2 has been included in the leads (indicated with dashed lines in the sketches 
of the cavities in the figure), in order to cause direct reflection. The magnetic 
field was increased as much as to produce a cyclotron radius smaller than 
a typical size of the cavity, and about twice the width of the leads. The 
cyclotron orbits are drawn to scale on the left of Fig. [3]for low and high fields. 
In all cases, the agreement between the numerical solutions of the Schrodinger 
equation and our maximum-entropy model is, generally speaking, found to 
be good [TS]. We should remark that in Figs. 3(e) and 3(f) four subintervals 
(treated independently) of 50 energies each had to be used, as each subinterval 
showed a slightly different optical S. This observation will be relevant to 
the discussion to be presented below on the validity and applicability of the 
information-theoretic model. 
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In order to elaborate further on the last remark, we investigated more 
closely a number of systems where the conditions for the approximate va- 
lidity of stationarity and ergodicity are not fulfilled. Fig. IHa) shows the 
statistical distribution of the conductance for a structure consisting of half 
a Bunimovich stadium coupled to a larger stadium. The system is time- 

600 



400 
200 





0.2 0.4 0.6 0.8 1 

T 

Fig. 4. (a) The conductance distribution obtained from the numerical integration 
of the Schrodinger equation for the structure shown in the inset and described in 
the text, compared with the prediction of Poisson's kernel, Eq. (|17p . with (S) 
extracted from the numerical data. The system is time-reversal invariant, so that 
we are dealing with the universality class P — 1. The agreement is reasonable, 
although one observes systematic deviations, as evidenced by the statistical error 
bars, (b) The square of the scattering wave function for a fixed energy inside the 
system considered in the analysis. We interpret the concentration of the wave 
function along the wall of the small cavity as a "whispering gallery mode" . (From 
Ref. [20].) 

reversal invariant, so that we are dealing with the universality class P — 1. 
The smaller stadium supports a "whispering gallery mode" -providing the 
short path in this problem- as illustrated in Fig. IHb), which shows the square 
of the absolute value of the scattering wave function for a fixed energy. The 
larger stadium provides a "sea" of fine structure resonances and should be 
responsible for the longer time scale in the problem. The two waveguides 
support = 1 open channel each and the system is time-reversal invariant 
(/? = !) (Ref. [20]). Statistics were collected by: 

i) sampling in an energy window AE. It was found that even for AE^s con- 
taining only sa 20 points -for which the S'(i?)'s were approximately uncorrelated- 
the energy variation of {S{E)) inside AE was not negligible. 

ii) constructing an ensemble of 200 structures by moving the position of 
the circular obstacle (shown in Fig. [4]) located inside the larger stadium. 

Again, the value of the optical S was extracted from the data and used 
as an input for Poisson's kernel, Eq. P7|) . We observe that the agreement 
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between theory and the computer simulation is reasonable. However, judging 
from the statistical error bars (and the results for other intervals as well), the 
discrepancies shown in the figure seem to be systematic. We interpret the 
situation described in i) above as giving evidence of two not very widely 
separated time scales in the problem. 

With the idea of decreasing the energy variation of {S{E)) inside the 
energy interval AE used to construct the histogram, AE was divided into 
two equal subintervals, giving the two conductance distributions shown in 
Figs. [SJa) and (b). We notice the improvement in the agreement between 




Fig. 5. The conductance distribution for the same structure as in Fig. U] but 
using energy intervals twice as small for the construction of the histograms. Panels 
(a) and (b) show the same data of Fig. [l] but analyzed inside each of the two 
subintervals. Panels (c) and (d) show the data for two other similar subintervals. 
The agreement with theory has improved considerably. (From Ref. |20j.) 

theory and numerics. The results for two other similar subintervals are shown 
in Fig. [HKc) and (d), where a similar conclusion applies. We remark that, in 
contrast to the results shown in Fig. [31 now we are "mainly" averaging over 
the ensemble of obstacles described in ii) above. 

We wish to remark that, no matter how much AE needs to be decreased 
in order to reduce the energy variation of {S{E)) across it, the improvement 
between theory and experiment would not be surprising if the number of 
resonances inside AE were kept constant, so as to keep approximate station- 
arity. This could be achieved, for instance, by increasing the size of the big 
cavity. However, this is not what has been done: the size of the big cavity 



QUANTUM CHAOTIC SCATTERING 13 



has been kept fixed, so that decreasing AE the number of resonances inside 
it decreases, and stationarity and ergodicity are ever less fulfilled. Yet, the 
agreement improves considerably: this we find surprising. 

The results of a recent calculation [21_ are even more intriguing. The 
conductance distribution was obtained numerically for the same structure as 
in the previous two figures, but now reducing the energy window AE literally 
to a point -so that the distribution actually represents statistics collected 
across the ensemble at a fixed energy- and the histogram was compared with 
the theoretical model. The results are shown in Fig. [S] for two different 
energies: they show a very good agreement with theory. 



100 




T T 



Fig. 6. The conductance distribution for the same structure as in Figs. |3]and[Sl but 
using an energy window AE = for the construction of the numerical histogram, 
which then represents statistics collected across the ensemble of obstacles at a fixed 
energy. Panels (a) and (b) correspond to two different energies. The agreement 
with theory is, in each case, very good. 

In Fig. [7| we present the conductance distribution for the same two- 
stadium structure as in the previous figures, again for the universality class 
P — 1, but now at an energy such that the waveguides support N — 2 open 
channels. Just as in Fig. [6l the numerical analysis was done collecting data 
across the ensemble for a fixed energy. The agreement, although not perfect, 
is reasonable. 

6.1 Discussion 

The theoretical model that we have used to compare with computer simu- 
lations, i.e., Poisson's kernel of Eq. (fT7|) . is based on the properties of an- 
alyticity, stationarity and ergodicity, plus a maximum-entropy ansatz. The 
analytical properties of the S matrix are of general validity. On the other 
hand, ergodicity is defined for stationary random processes, and stationarity, 
in turn, is an extreme idealization that regards S{E) as a stationary random 
(matrix) function of energy: stationarity implies, for instance, that the op- 
tical {S{E)) is constant with energy and the characteristic time associated 
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T 

Fig. 7. The conductance distribution for the structure shown in the inset, at an 
energy such that N = 2 channels are open. The system is time-reversal invariant, 
so that we are dealing with the universality class /3 = 1. The energy interval AE 
was literally reduced to a point (after Ref. [20]). The agreement between theory 
and the simulation is reasonable. 

with direct processes is literally zero. Needless to say, in realistic dynamical 
problems like the ones considered in the present section, stationarity is only 
approximately fulfilled, so that one has to work with energy intervals AE 
across which the "local" optical {S{E)) is approximately constant, while, at 
the same time, such intervals should contain many fine-structure resonances. 
Above we saw examples in which this compromise is approximately fulfilled, 
and cases in which it is not. Yet, we have given evidence that reducing the 
size of AE, even at the expense of decreasing the number of resonances inside 
it, improves significantly the agreement between theory and numerical exper- 
iments. We went to the extreme of reducing AE literally to a point, which 
amounted to leaving the energy fixed and collecting data over the ensemble 
constructed by changing the position of the obstacle: in the example shown 
in Fig. [n]we found an excellent agreement of the resulting distribution with 
Poisson's kernel. 

An essential ingredient to construct Poisson's kernel of Eq. (|17p is the 
reproducing property of Eq. (|16p . This is a property of the ensemble, defined 
for a fixed energy: at present we only know how to justify it through the 
properties of stationarity and ergodicity. However, the results we have shown 
(Fig. [5l which is "almost" an ensemble average, and Fig. [5] which is literally 
an ensemble average), seem to suggest that Poisson's kernel is valid beyond 
the situation where it was originally derived: i.e., it is as if the reproducing 
property of Eq. ()16|) were valid even in the absence of stationarity and 
ergodicity. We do not have an explanation of this fact at the present moment. 

7 Conclusions 

In this contribution we reviewed a model that was originally developed in 
the realm of nuclear physics, and applied it to the description of statistical 
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Quantum Mechanical scattering of one-particle systems inside cavities whose 
underlying classical dynamics is chaotic. The same model is also applicable 
to the study of scattering of classical waves by cavities of a similar shape. 

In the scheme that we have described, the S{E) matrix is idealized as a 
stationary random (matrix) function of E, satisfying the condition of ergod- 
icity. It is assumed that we have two distinct time scales in the problem, 
arising from a prompt and an equilibrated component: the former one is as- 
sociated with direct processes, like short paths, and is described through the 
optical matrix (S). 

The statistical distribution of S is proposed as one of maximum entropy, 
once the reproducing property described in the text is fulfilled and the system- 
dependent optical (5) is specified. In other words, system-specific details 
other than the optical matrix (S) are assumed to be irrelevant. 

The predictions of the model for the conductance distribution in the pres- 
ence of direct processes should be experimentally observable in microstruc- 
tures where phase breaking is small enough, the sampling being performed 
by varying the energy or shape of the structure with an external gate voltage. 
They have also been observed in microwave scattering experiments. 

Comparison of the theoretical predictions with computer simulations indi- 
cate that the present formulation has a remarkable predictive power; indeed, 
it describes statistical properties of the quantum conductance of cavities, like 
its average, its fluctuations, and its full distribution in a variety of cases, both 
in the absence and in the presence of a prompt response. Here we have con- 
centrated on the conductance distribution, mainly in the presence of prompt 
processes. 

In Ref. |15j we found that, with a few exceptions, it was possible to 
find energy intervals AE containing many fine-structure resonances in which 
approximate stationarity could be defined; the agreement was generally found 
to be good. In contrast, in Ref. [2D] we encountered cases in which this was 
not possible and we undertook a closer analysis of the problem. We found that 
reducing AE, even at the expense of decreasing the number of fine-structure 
resonances inside it, improved the agreement considerably. In a number of 
cases, AE was reduced literally to a point and the data were collected over 
an ensemble constructed by changing the position of an obstacle inside the 
cavity: in other words, Poisson's kernel was found to give a good description 
of the statistics of the data taken across the ensemble for a fixed energy. As 
a result, Poisson's kernel seems to be valid beyond the situation where it was 
originally derived, where the properties of stationarity and ergodicity played 
an important role. It is as though the reproducing property of Eq. p^ . 
which is a property of the ensemble defined for a fixed energy, were valid 
even in the absence of stationarity and ergodicity, which were originally used 
to derive it. At the present moment we are unable to give an explanation of 
this fact. 
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